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Motivated by recent lattice QCD studies, we explore the effects of interactions on strangeness 
fluctuations in strongly interacting matter at finite temperature. We focus on S-wave Ktt scattering 
and discuss the role of the ifo(800) and K*(1430) resonances within the S-matrix formulation 
of thermodynamics. Using the empirical Ktt phase shifts as input, we find that the Ktt S-wave 
interactions provide part of the missing contribution to the strangeness susceptibility. Moreover, it 
is shown that the simplified treatment of the interactions in this channel, employed in the hadron 
resonance gas approach, leads to a systematic overestimate of the strangeness fluctuations. 

PACS numbers: 24.10.Pa, 25.70.Bc, 25.70.Ef, 25.75.-q 


I. MISSING INTERACTION STRENGTH IN 
STRANGE SECTOR 

A recent study, comparing QCD thermodynamics ob¬ 
tained on the lattice with the hadron resonance gas 
(HRG) model [T], indicates that additional interaction 
strength, beyond that embodied by well established 
strange resonances [2], may be needed in the HRG 
model to remove disparities with the lattice results. 
In particular, the HRG results for the strangeness and 
mixed strangeness-baryon number susceptibilities (xss 
and xbs) are clearly below those of the lattice, while the 
results for the thermodynamic pressure and the baryon 
number susceptibility are in good agreement. 

This motivates the search for hitherto unknown 
strange hadrons, which could reduce or eliminate this dis¬ 
crepancy. In the PDG database, there are around twenty 
unconfirmed states with a mass below 2.0 GeV. Although 
these are not established resonances, the interactions in 
the corresponding scattering channels may yield impor¬ 
tant contribution to thermodynamic quantities. 

More generally, a possible origin of the discrepancy is 
interaction strength in channels carrying net strangeness 
that so far have not been accounted for. Given the 
corresponding empirical scattering phase shifts, both 
confirmed and unconfirmed resonances as well as non¬ 
resonant interactions can be handled in a unified, model- 
independent way, using the S-matrix approach of Ref. [3] . 

The strange scalar channel, with the unconfirmed 
Arg(800) resonance, a.k.a. k, is a prime candidate. Since 
the corresponding phase shifts for S-wave Ktt scattering 
are fairly well determined, this channel is well suited for 
the S-matrix approach. In addition, the counterpart of 
K in the scalar-isoscalar channel, the /o(500), a.k.a. a, 
though considered to be established [4] , is unlike a typical 
resonance. Since the tttt S-wave phase shifts are known 
with reasonable accuracy, also this channel is a prime 
candidate for the S-matrix approach to thermodynam¬ 
ics. In this study we focus on the strange scalar channel 
and its contribution to strangeness susceptibilities. 


With the relatively low mass of the interaction strength 
in the k channel, it potentially has a large impact on the 
thermodynamics, in particular on xss, owing to the mod¬ 
erate suppression by the Boltzmann factor. In Fig. 
we illustrate the effect of the k resonance on pressure 
and strangeness fluctuation within the HRG approach. 
For the PDG particle spectrum, we use only confirmed 
baryons (i.e. three and four star resonances) and es¬ 
tablished mesons. The contribution of k to the ther¬ 
modynamics is approximated by that of an ideal gas of 
zero-width mesons with mass = 0.682 GeV and de¬ 
generacy four. Indeed, the inclusion of this single state 
improves the HRG result on xss dramatically, while the 
agreement in the thermodynamic pressure persists. How¬ 
ever, owing to the fairly large width, the treatment of 
the K resonance as a zero-width particle is questionable. 
Consequently, a systematic approach, where all interac¬ 
tion effects are treated consistently, is called for. 

In this paper, we assess the effect of interactions in the 
K channel on the thermodynamics using the S-matrix ap¬ 
proach [5] . For elastic scattering, the resulting expression 
reduces to the Beth-Uhlenbeck form for the second virial 
coefficient, expressed in terms of the scattering phase 
shift |S]. In Ref. [5] this scheme was applied to com¬ 
pute the contribution of ttN interactions to the baryon 
number and the tt transverse momentum spectrum in 
hadronic matter at moderate temperatures and densities. 
The method yields an effective spectral weight, which is 
relevant for the partition sum and thus allows one to 
compute the interaction contribution to various thermo¬ 
dynamic observables. 

The paper is organized as follows. In section H we 
describe the parametrization of the empirical Ktt phase 
shifts. In particular, we implement the constraint pro¬ 
vided by the empirical scattering length in this channel. 
Moreover, we discuss the connection between the phase 
shifts and thermodynamic quantities. In section HI we 
apply this formalism to study the effect of interactions in 
the K channel. We compare these results with those of the 
standard HRG approach, where the k meson is treated 
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FIG. 1: (color online). Left: The thermodynamic pressure (normalized to T'^) computed in HRG using only established 
resonances (PDG, broken dashed line) and adding the unconfirmed k (PDG + kappa, full line). The lattice results on pressure 
are from Refs. [7l[8]. The dashed line shows the results where the k channel is treated within the S-matrix formulation (see 
text). Right: The corresponding results for strange susceptibility xss (normalized to T^). The lattice results on fluctuations 
are from Ref. [9] [T0| . 


as a Breit-Wigner resonance with an energy-independent 
width and study the influence of width and scattering 
length on strangeness fluctuations. In the final section, 
we present our conclusions. 


clear from the above expression that the self-energy will 
develop an imaginary part when the invariant four mo¬ 
mentum exceeds the threshold, i.e. s = k'^ > 4to^. The 
width of $ in this model reads 


II. S-WAVE SCATTERING AND THE 
S-MATRIX APPROACH 

The scattering phase shift contains the necessary phys¬ 
ical information to study resonances using scattering 
data. Since the k meson has the quantum numbers 
I{J^) = 1 ( 0 ''"), the relevant phase shift to consider 
is that of kaon-pion scattering in the S-wave, isospin 
1 = 1/2 channel (Jq' ). We shall begin by collecting some 
basic field theoretical results to establish the connection 
between resonance width and the scattering phase shift 
pertinent to the study of the k resonance. 


A. S-wave decay of resonance 

Consider the decay of a scalar particle $ —>■ (j) + (j), 
through an interaction term Ci = —g . The self¬ 
energy of $ is, to leading order in perturbation theory, 
given by 




2ig^ 


dH 1 1 

(27r)4 l'^ — m? {I — fc)2 — 'm? 



a:(l — x) k'^) tt], 

( 1 ) 


where m is the mass of the particle 4>. The decay rate 
is obtained from the imaginary part of E$(fc^). It is 


VS 


= 0[s - (2m)2] Vl-(2m)Vs^. (2) 


In the more general case of S-wave decay of a resonance 
into two particles with different masses mi and m 2 , one 
finds 


/ N « 21 Pcm{s) 

7(s) = w ^ s - mth] - 

2 s 

PcAiis) = m\Js\Jl- (3) 

Am = mi — m 2 
mth = mi -I- m2. 


Here we have introduced the notations a = g^/dTr and the 
center of mass momentum Pcm{s). Note the different 
symmetry factor in Eq. (§, owing to the distinguishable 
particles in the final state. As we show below, the energy 
dependence of the width 7 ( 5 ) is crucial for reproducing 
the S-wave phase shifts near threshold. 


B. Parametrization of the S-wave Kn phase shift 

Although the expression for the decay width in Eq. (§ 
is obtained from a perturbative one-loop calculation, it 
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a^GeV") 

M«(GeV) 

aK„-(GeV^) 

Mk* (GeV) 

r^F2(GeV-i) 

3.0098 

0.905 

1.437 

1.41 

3.57 


0.81 


TABLE I: Parameters used to model the Ktt scattering phase 
shifts in the S-wave, isospin 7 = 1/2 and 1 = J> 12 channels 

[HI- 



1=1/2 

CLq UT-tt 

1=3/2 
tto m^r 

Ishida et al. [TT] 

0.393 

-0.112 

Biittiker et al. [13] 

0.224(22) 

-0.045(8) 


TABLE II: The S-wave K-k scattering lengths obtained with 
the model m confronted with an empirical value. 

provides a general form for parametrizing the phase shifts 
of S-wave scattering. We account for the contribution of 
the two lightest 0+ strange resonances, k and 7fQ(I430), 
to the Ktt phase shifts by using an energy-dependent 
Breit-Wigner form 

,5re8(s) = tan~^( 

= (4) 

for each resonance. Here a and Mq are free parameters, 
which are fitted to the data. 

However, in addition to the resonance contribution, a 
repulsive background contribution is needed for a suc¬ 
cessful description of the empirical phase shifts in this 
channel. Following m, we parametrize the background 
using the phase shift of a hard sphere 

Sbg{s) = —rcPcAiis), (5) 

where Tc is the radius of the repulsive core. The total 
phase shift is given by the sum of the resonance 
contributions and the background 

+ Sk» + Sbg- ( 6 ) 

Using the parameters obtained by Ishida et al. [n (sum- 
marized in Table [I]), Eqs. Q - (|^ provide a good descrip¬ 
tion of the experimental data up to 1.6 GeV, as shown in 
Fig.H 

1 /2 

An additional constraint on the fit of 5q comes from 
the scattering length [14], which is related to the 
phase shift near the threshold by 

^ mth) := aJ^PcMis) + ( 7 ) 




FIG. 2: (color online). Top: The Ktt scattering phase shift in 
the S-wave, isospin 7 = 1/2 channel. The experimental results 
are obtained from Ref. [nisi. The solid line corresponds 
to the parametrization discussed in Eqs. @ - § using the 
parameters depicted in Table [I| The band corresponds to dif¬ 
ferent values of the scattering length, obtained by adjusting 
the regulator in Eqs. ([^ The dashed line corresponds to 
the phase shift with = O.lSmj^^. Bottom: Similarly for 
the S-wave, isospin 7 = 3/2 channel. In this case, the dashed 
line corresponds to the phase shift with = 0.045 


However, the threshold behavior is not uniquely deter¬ 
mined by the data. Hence, additional input is needed to 
obtain an accurate description of the scattering length. 

The value for the 7=1/2 S-wave Ktt scattering length 
obtained in the model is high compared to that obtained 
in a dispersive analysis oi Ktt scattering [TS] (see Table 
H). Low-energy theorems based on the current algebra 
and the partially conserved axial-vector current (PCAC) 
predict a lower value of « 0.14 to”^ At 

next-to-leading order in chiral perturbation theory [18) . 
the scattering length is « O.ISto”^, while agreement with 
the dispersive approach can be obtained at NNLO [19]. 
The extraction of the Ktt scattering length from lattice 
QCD is at present not conclusive. So far, such calcula¬ 
tions were done with large pion masses, resulting in large 
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values of the scattering lengths EHIIII]. The extrapola¬ 
tion to physical value of the pion mass is delicate. Using 
chiral perturbation theory at next-to-leading order for 
the extrapolation, Fu [H] finds a somewhat low value, 

1 /2 

To cover the range of uncertainty in , we introduce 
a regulator 


In this approach, the thermodynamic potential of an 
interacting system of pions, kaons and resonances is, to 
leading order, given by the sum: 

n = rijr + + ^int- (H) 

The first two terms are the ideal gas expressions for pions 
and kaons: 


J-(s) = 1 - 


/o 

1 -b s/A2 


( 8 ) 


in the K-contribution to the phase shift such that 


(5„(s) = tan ^ 


J-(s) 


-v/s7(s) ' 
s - M2 _ ■ 


(9) 


Numerically we use A = 0.381 GeV and vary /o to ob¬ 
tain a scattering length between 0.18 for /o = 1 and 
0.4 for /q = 0. 

Before we end the discussion of phase shift, we com¬ 
ment on two important features. First, the shape of 
( 5 q' (s) differs qualitatively from that of a narrow res¬ 
onance. In the limit of vanishing width, as usually as¬ 
sumed in the HRG model, the phase shift would become 
a step function which reaches the value of 180° at a mass 
of y/s « 0.682GeV. A comparison with Fig. [^clearly 
shows that the k meson cannot be treated as a narrow 
resonance. Second, the behavior of the phase shift at 
threshold is determined by the orbital angular momen¬ 
tum. In an S-wave, the derivative of the phase shift with 
respect to s 


= STU J ^ {ln[l - (^3) 

nK = 2 TV [ (ln[l-e-'^(Vp^+^+^s)] 

-b ln[l - I ^ ( 13 ) 

where /ig is the strangeness chemical potential and the 
two terms in Hk are due to kaons and antikaons, respec¬ 
tively. Finally, the last term in Eq. © accounts for Ktt 
interactions. In the HRG approach, Hint is given by the 
sum of all relevant resonances treated as an ideal gas of 
stable particles: 


= E 2?^^/ ^ {ln[l - 

-b ln[l - I (14) 

The degeneracy factor accounts for the two possible 
isospin states in the I{J^) = 5 ( 6 “'') channel. The ther¬ 
modynamic pressure is computed by 


— ~ Pcm{s ( 10 ) 

diverges at threshold, due to the fact that Pqm = 
dPcm/ds diverges at s = (see Eq. @)- As we dis¬ 
cuss in the next section, this fact determines the behav¬ 
ior of the effective spectral weight derived from the phase 
shifts. 


C. S-matrix approach and thermodynamics 

Given the parametrization of the Ktt S-wave phase 
shifts presented above, we are now ready to formulate 
the thermodynamics. The tool of choice is the S-matrix 
approach , which provides a systematic way to account 
for interactions in a many-body system in thermal equi¬ 
librium. The leading order correction, which is deter¬ 
mined by the two-body scattering phase shift, is equiv¬ 
alent to the second virial coefficient [5]. We apply this 
formalism to compute the interaction contribution of Ktt 
scattering in the k, channel to the thermodynamics of 
strongly interacting matter in the hadronic phase. 


P = - 


Vt 

V' 


(15) 


Another key quantity of interest is the strangeness sus¬ 
ceptibility, which is obtained by taking derivatives of the 
thermodynamic pressure with respect to the strangeness 
chemical potential 


d^P 

XSS = ^— 

o^^so^^s 

The results for the thermodynamic observables in HRG 
are shown in Fig. 

In the S-matrix approach, the interaction contribution 
to the thermodynamic potential involves an integral over 
the invariant mass M = 


(16) 



X {ln[l - (17) 

-b ln[l - . 



















5 




FIG. 3: (color online). Left: Comparison of the weight function B{M) (solid line) and the double Breit-Wigner spectral function 
for K and 7Fo(1430) (dashed line), both in units of GeV“^. The band corresponds to different values of the scattering length, 
bounded by the solid and dashed lines, which corresponds to the scattering length of = 0.4 and 0.18 respectively. 
Right: The interaction contributions to the strangeness susceptibility in the S-wave 7=1/2 channel, obtained using different 
spectral weights. 


with the effective weight function I3E] 


III. INFLUENCE OF THE k CHANNEL 


B{M) = 2 


d 

dM 


6{M), 


which satisfies the normalization condition 


(18) 



dM 


B{M) = 1, 


(19) 


provided the phase shift has the property S{oo) —>■ tt. 
We note that the weight function B is in principle de¬ 
fined in any channel, irrespective of the existence of a 
corresponding resonance. Moreover, even for a well de¬ 
fined resonance, the weight function differs from the cor¬ 
responding spectral function |5]. The two functions are 
identical only in the limit, where the width of the reso¬ 
nance vanishes. 

Using the thermodynamic potential 0 , one can com¬ 
pute the interaction contribution to thermodynamic ob¬ 
servables. Thus, e.g. the interaction contribution to the 
thermodynamic pressure is given by 


/■°° dM 

P^ut = / —B{M)Pt{M), (20) 

Jmth 

where 

PriM) = -2j ^ |ln[l - ( 21 ) 

-b ln[l-e-v'p'+^'+'^®]| 

with P = PfT'^, p = p/T, M = M/T and fis = ps/T. 
The interaction effects on the strangeness susceptibility 
will be discussed in the following section. 


The parametrization of the phase shifts and weight 
function presented above allow for an assessment of var¬ 
ious approximate descriptions of the interaction effects 
in the Ktt channel. For example, a standard Breit- 
Wigner resonance with an energy-independent width 
7(5) —>■ 7bw is easily accommodated by neglecting the 
s-dependence of the numerator ■Psj(s) in the phase shift 
formula Eq. @ 0 , thus 


W = ^^S{M) 

2M 


2M7BW 


(M2 - M2)2 + M^-ff 


( 22 ) 


BW 


The weight functions in different approximation 
schemes are shown in Fig. (left panel). The valid¬ 
ity of the weight function is limited to energies below 
M = 1.6 GeV, the highest energy included in the fit of 
the phase shifts. For the computation of thermodynamic 
observables, e.g. pressure in Eqs. (pOl) - (pTl), the integral 


over M converges well below this energy for tempera¬ 
tures up to 0.16 GeV, owing to the suppression by the 
Boltzmann factor. 

One characteristic feature of the weight function B for 
an S-wave channel is that it diverges at the threshold, 
as seen in the left panel of Fig. This singularity is, 
however, integrable, and its sign as well as its strength 
are directly related to the scattering length in the cor¬ 
responding partial wave. A lower value of the scattering 
length tends to reduce the strength of the weight function 
near the threshold, shown as the blue band of Fig. The 
standard Breit-Wigner form, on the other hand, does not 
exhibit such a divergence. 
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FIG. 4: (color online). Left: Comparison of the S-wave weight functions (M) (GeV“^) in different isospin channels and their 
weighed sum. Right: The interaction contributions to the strangeness susceptibility (normalized to T^) from the I = 1/2 and 
1 = 3/2 S-wave channels. 


In addition, in Fig. (left panel), we observe a sec¬ 
ondary peak appeared near 1.4 GeV. This clearly corre¬ 
sponds to the iLg(1430) resonance. 

In the right panel of Fig. we show the dependence 
of XSS on the weight function. We observe that the S- 
matrix approach yields a result that lies between those 
obtained within the standard Breit-Wigner, for Ffg (1430) 
with and without k. Moreover, the dependence of xss on 
the scattering length is displayed. We see that a larger 
positive scattering length provides more support in the 
low mass region and hence gives an enhanced contribu¬ 
tion to the susceptibility. 

We now discuss the origin of the suppression of thermo¬ 
dynamic observables when treating the interaction based 
on K-channel phase shift. Previous studies have stressed 
the importance of using B instead of the standard spec¬ 
tral function [5]. For the case of A resonance, where 
P-wave scattering is involved, the B function tends to 
enhance the low mass contribution to the thermodynam¬ 
ics and results in an overall increase in the observables 
beyond those treated by the standard Breit-Wigner ap¬ 
proach. In the current study involving an S-wave scatter¬ 
ing, the enhancement effect near the threshold is, how¬ 
ever, compensated by the relatively slow increase in the 
phase shift below 1.3 GeV. Unlike a typical resonance, 
the phase shift in the k channel does not reach 180° 
before iFo(1430) emerges. The slow rise of the phase 
shift in the low mass region limits the strength of weight 
function B. Consequently, the thermodynamic observ¬ 
ables calculated in the S-matrix approach for k becomes 
strongly suppressed, making its contribution too small 
to remove the disparity between HRG and lattice results 
(see Fig. [^. 

Thus, our calculation shows that the contribution of 
Ktt interactions to the strange susceptibility is substan¬ 
tially lower in the current consistent treatment compared 
to the HRG description. As a result, the inclusion of 


the K channel does not resolve the issue of the missing 
strength in the strange sector. Clearly, a careful analysis 
of the interaction strength in other strangeness carrying 
channels is called for. 


IV. THE EFFECT OF J = 3/2 Kn SCATTERING 

For completeness, we also assess the contribution of 
the isospin I = 3/2 S-wave channel to the strangeness 
fluctuations. As pointed out in Ref. [52], the inclusion of 
this channel partly cancels the effect of the isospin I = 
1/2 channel. This is expected since the phase shift in the 
I = ^12 channel, and consequently the weight function 
H^=3/2(M), are negative, corresponding to a repulsive 
interaction. 

Given that this channel involves only non-resonant Ktt 
scattering, we employ the repulsive core expression for 
the phase shift 


= -ri=^^^PcM{M). 


(23) 


The fit parameter ri = 0.112suggested by 
Ishida et al. El, yields a high value for the scatter¬ 
ing length compared to that obtained in the dispersive 
analysis of m (see Table II). To cover the range of uncer¬ 
tainties in the scattering lengths, we employ the regulator 
introduced in Eq. ([^. We use A = 0.381 GeV and dial 
/o to obtain a scattering length between 0.045 for 
/o = 2.28 and 0.112 to"! for /o = 0. 

The resulting phase shift 5^^ is shown in Fig. (lower 
panel). It is evident that Eq. (23) cannot capture the 
full features of the phase shift, and a more refined ap¬ 
proach (like e.g. [23]) can achieve a more satisfactory 
description of the non-resonant scattering. Nevertheless, 
the parametrization (23) is sufficient for our current dis¬ 


cussion. 
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The weight function S, in the thermodynamic poten¬ 
tial discussed in Eqs. 0 - pTj ), is then modified as 
follows: 


(24) 

where (M) is the weight function in isospin channel I 
and the factor of 2 in front of accounts for the 

relative isospin degeneracy factor of the I = 2> 12 and 
7 = 1/2 channels. 

The corresponding contributions to the weight func¬ 
tion and the strangeness susceptibility are shown in 
Fig-i The partial cancellation between the two isospin 
channels is evident. As a result, the enhancement of 
strangeness fluctuations due to S-wave K-k scattering is 
reduced by 70%. This effect, which is not accounted for 
in the HRG model, lends further support to our conclu¬ 
sion that a consistent treatment of low-mass resonances 
requires a careful analysis, including also non-resonant 
interactions. A natural framework for such studies is of¬ 
fered by the S-matrix approach employed in this paper. 


V. CONCLUSION 

This study set out to explore possible sources of miss¬ 
ing strength in the strangeness susceptibility, suggested 
by lattice results. The 77 q( 800) resonance, a.k.a. k, 
which is not an established resonance in the PDG com¬ 
pilation, appears to be a promising candidate. Indeed, 


within the treatment of the hadron resonance gas (HRG) 
model, we found that this single state alone accounts for 
the missing contribution in the strange susceptibility. 

However, owing to the large width of the n meson 
and the significant non-resonant background, the HRG 
model does not provide an accurate description of the 
interaction contributions. In fact, a consistent treatment 
of all Ktt S-wave interactions within the S-matrix ap¬ 
proach shows that a simplihed (HRG) treatment of the 
interactions in these channels, using a Breit-Wigner spec¬ 
tral function for each resonance and ignoring the non¬ 
resonant background, systematically overestimates the 
contribution to strangeness fluctuations. 

In summary, the Ktt S-wave interactions provide only 
a part of the missing contribution to the strangeness 
susceptibility, indicated by recent lattice QCD results. 
Whether the remaining discrepancy can be resolved by a 
consistent treatment of other strangeness carrying chan¬ 
nels will be explored in future studies. 
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